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ABSTRACT 

MSH 11-62 (G291T— 0.9) is a composite supernova remnant for which radio and X-ray observations have identified 
the remnant shell as well as its central pulsar wind nebula. The observations suggest a relatively young system 
] expanding into a low density region. Here we present a study of MSH 11-6^ using observations with the Chandra, 
XMM, and Fermi observatories, along with radio observations from the Australia Telescope Compact Array (ATCA). 
We identify a compact X-ray source that appears to be the putative pulsar that powers the nebula, and show that the 
] X-ray spectrum of the nebula bears the signature of synchrotron losses as particles diffuse into the outer nebula. Using 
data from the Fermi LAT, we identify 7-ray emission originating from MSH 11-62. With density constraints from 
the new X-ray measurements of the remnant, we model the evolution of the composite system in order to constrain 
the properties of the underlying pulsar and the origin of the 7-ray emission. 

Subject headings: radiation mechanisms: non-thermal — ISM: supernova remnants: — ISM : indi- 
vidual (MSH 11-6^) — stars: neutron 



1. INTRODUCTION 

The creation of an energetic pulsar in a supernova ex- 
plosion often results in the formation of a composite su- 
pernova remnant (SNR) that can be characterized by 
emission from the remnant, an associated pulsar wind 
nebula (PWN), and the neutron star (NS) itself (see, e.g., 
Gaensler & Slane 2006). Recent studies of such systems 
have resulted in discoveries of new pulsars, constraints 
on the spectra of particles injected from the pulsars into 
their nebulae, and evidence for interactions between the 
SNR reverse shock (RS) and the PWN. Broadband ob- 
servations of composite SNRs have been crucial for im- 
proving our understanding of the detailed evolution of 
I these systems, particularly with the vast improvements 
in capabilities provided by current X-ray and 7-ray ob- 
servatories. 

; MSH 11-6^ (Mills, Slee, & HiU 1961) is a centrally- 
brightened SNR whose properties strongly suggest the 
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presence of a central nebula powered by an active pul- 
sar. Radio observations reveal a remnant with a dis- 
tinct bar-like central structure whose flat radio spectrum 
{ar — 0.29 ± 0.05, where Si, oc j/~"' ) is similar to that 
observed for PWNe, and polarization measurements at 5 
and 8.4 GHz indicate that the magnetic field in the cen- 
tral region of the SNR is aligned with the long axis of the 
bar (Roger et al. 1986); the flux density at 8.4 GHz is 
10.4 ± 0.4 Jy. No radio pulsar associated with MSH 11- 
62 has been detected. An 11 hr observation at 1.4 GHz 
failed to detect any pulsations from this region, down 
to a sensitivity limit of 0.028 mJy assuming a 10% duty 
cyclO 

The distance to MSH 11-62 is poorly known. Limits 
on HI absorption indicate a distance greater than 3.5 kpc 
(Moffett, Gaensler, & Green 2001; Moffett et al. 2002). 
We adopt a fiducial distance of 5 kpc in derived quantities 
below, and scale values to d^, the distance in units of 
5 kpc. Based on the radio synchrotron spectrum, the 
minimum-energy magnetic field in the central nebula has 

been estimated as i? ^ 65rf5 ^^^(0p/2'.5)~^/^/iG, where 
6'p is the PWN radius in arcmin (Harms, Hughes, & 
Slane 1998). 

X-ray emission from MSH 11-62 was first identified by 
Wilson (1986) using the Einstein Observatory. X-ray ob- 
servations with the ASCA observatory reveal a spectrum 
best described by a combination of thermal and non- 
thermal models (Harms, Hughes, & Slane 1998). The 
hard X-ray emission is spatially concentrated in the cen- 
tral regions, providing further evidence for the presence 
of a PWN within the remnant, but the limited angular 
resolution of ASCA precluded detection of the pulsar it- 
self. The derived spectral index for the X-ray power law 
component is F = 2.0to'3 ^^'^ the associated luminosity 
is L^(0.2 - 4 keV) ~ 6 x lO^^dj erg s'^. The best-fit 
temperature of the thermal component, presumably as- 
sociated with emission from shock-heated gas in the SNR 

^■^ F. Camilo, private communication. 
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Fig. 1.— ATCA 1.4 GHz image of MSH 11-62. The remnant is 
characterized by a faint shell with a bright central bar. Contours 
levels are 3.6 and 1.5 mjy beam^^. The cross marks the position 
of the X-ray point source described in Section 3. 

shell, is kT ^ 0.8 keV. These observations clearly imply 
that MSH 11-6^ harbors a young pulsar that powers a 
central PWN. Chandra observations, which we discuss 
in detail below, reveal a compact X-ray source that is 
almost certainly the pulsar counterpart (Harms et al. 
2002; Kargahsev & Pavlov 2008). 

Observations with the Compton Gamma-Ray Observa- 
tory discovered the source 3EG J1102— 6103, for which 
MSH 11-6^ falls within the position error circle. The flux 
is F(> 100 MeV) = (3.3±0.6) x lO"'' photons cm"^ s"!, 
with a power law index F = 2.5 ± 0.2 (Hartman et al. 
1999). This source was suggested as a 7-ray counterpart 
to MSH 11-6^ (Sturner & Dermer 1995), but was also 
subsequently suggested as the counterpart to the nearby 
young pulsar J1105— 6107 (Kaspi et al. 1997). Given the 
large uncertainty in the position of the 7-ray source, the 
origin of its emission remained uncertain. 

Here we report on observations with the ATCA, the 
Chandra X-ray Observatory, the XMM-Newton X-ray 
Observatory, and the Fermi Camma-Ray Space Telescope 
that provide sufficient sensitivity and resolution to iden- 
tify the neutron star in MSH 11-6^, and to provide a 
much-improved characterization of the properties of this 
composite system. In Sections 2 we discuss our ATCA 
observations. In Section 3 we summarize the Chandra 
and XMM observations, with detailed investigations of 
the PWN, neutron star, and SNR shell. In Section 4 we 
discuss observations with the Fermi LAT, and in Section 
5 we present discussions of the results along with our 
conclusions. 

2. RADIO OBSERVATIONS 

2.1. Observations and Data Reduction 

Radio observations of MSH 11-6^ were carried out 
with the ATCA, a 6-element interferometer located near 
Narrabri, New South Wales, Australia. Observations 
were carried out on three separate dates using three dif- 
ferent array configurations: on 1999 Feb 10 (750C array). 
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Fig. 2. — Chandra image of MSH 11—62. Colors correspond to 
different X-ray energy bands: red = 0.75 - 1.45 keV; green = 1.45 
- 2.58 keV; blue = 2.58 - 6.0 keV. The SNR shell is evident as 
soft (red) X-ray emission (the portion of the shell in the northeast 
falls outside the detector field of view). The central bar and asso- 
ciated compact source are white in color, corresponding to harder 
emission. The inset figure shows the central regions of the PWN 
containing the pulsar candidate (white) and the star Tycho 8959 
255 1 (red). 

1999 Feb 22 (6C array), and 1999 Mar 24 (1.5B array). 
In all cases, radio continuum data were recorded in all 
four polarizations, over a bandwidth of 128 MHz cen- 
tered at 1384 MHz. At each epoch, MSH 11-6^ was 
observed for 8-11 hours spread over a 12-hour synthesis, 
combined with regular short observations of the extra- 
galactic source PKS B1036~697 to calibrate atmospheric 
gain variations and polarization leakages. The bright 
source PKS B1934-638 was also observed at each epoch 
to provide an overall flux calibration. 

Following editing of the data and calibration in the 
MIRIAD package using standard techniques, the visibility 
data-set was imaged using multi-frequency synthesis and 
uniform weighting, and then deconvolved with a maxi- 
mum entropy algorithm. The final total intensity image 
has an angular resolution of 8'.'8 x 6'.'6, with the major 
axis oriented at a position angle 3° east of north, and 
an RMS sensitivity of 0.1-0.2 mJy beam^^. Preliminary 
investigations of these observations have been presented 
by Moffctt, Gaensler, & Green (2001) and Moffett et al. 
(2002). 

2.2. Radio Structure and Properties 

The ATCA 1.4 GHz image of MSH 11-6^ (Figure 
1) reveals a roughly circular morphology with a bright 
central bar-like structure that corresponds with the flat- 
spectrum region identified by Roger et al. (1986). Par- 
tial limb-brightening is observed along the northwestern 
and southeastern boundaries of the SNR shell. Consid- 
erable filamentary structures are observed both along 
the limb and in the SNR interior. The bright central 
bar reveals several regions of enhanced radio emission, 
with the brightest feature located near the southwest- 
ern end of the structure (at an approximate position of 
llhiim4gs _gQ03g/27//-) ^hc structure of the central bar 
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is similar to that observed in several other PWNe, in- 
cluding Gil. 2-0.3 (Tam et al. 2002), N157B (Lazen- 
dic et al. 2000), and G328.4+0.2 (Johnston, McClure- 
Griffiths, & Koribalski 2004; Gelfand et al. 2007). Its 
apparent similarity to the PWN in G328.4+0.2, in par- 
ticular, may suggest that the nebula has been compressed 
by the SNR reverse shock (RS), which may have prop- 
agated more quickly from the NW/SE directions, where 
the SNR limbs show slight brightness enhancements. 

The total 1.4 GHz flux from the SNR is 14 Jy, with 
4.8 Jy associated with the central nebula. This revised 
flux value results in a new minimum- energy estimate for 

the magnetic field: B « 36d^^/'^(6'p/2'.5)-^/VG. As 
we discuss in Section 5.3, however, some evolved PWNe 
are known to be particle-dominated, with magnetic field 
strengths far below the minimum-energy (or cquiparti- 
tion) values. 

3. X-RAY OBSERVATIONS 

A 50 ks Chandra ACIS observation of MSH 11-6^ was 
performed on 2002 April 8 (ObsID 2782). The central 
region of the SNR was placed on the S3 chip of the de- 
tector, and portions of the shell were imaged on the ad- 
jacent chip S2. The extent of the SNR (with a radius 
-RsNR ~ 5.7 arcmin) is considerably larger than the field 
covered by the S3 chip, which resulted in the northeast- 
ern regions of the remnant falling outside of the detector 
field of view. 

MSH 11-6^ was observed for 40 ks with XMM, on 
2002 February 6 (ObsID 0051550101). The M0S2 and 
pn detectors were operated in full frame mode with the 
MEDIUM filters, providing fuU coverage of the SNR, 
while MOSl was operated in timing mode with the THIN 
filter. 

3.1. Observations and Data Reduction 

The Chandra data were reduced using standard screen- 
ing prescriptions. All data were reduced using tools from 
CIAO 3.2.2.. After screening, the remaining exposure 
time was 49.4 ks. 

The Chandra observations show hard emission ex- 
tended along the radio bar (Figure 2), although the emis- 
sion does not extend as far to the NE as the radio bar. 
There is also a clear shell of soft X-ray emission that 
corresponds well with the morphology of the radio SNR, 
including a small intensity enhancement along the posi- 
tion of the western radio limb-brightening, and a large 
under luminous region to the northwest of the radio bar. 

Of particular significance in the Chandra image is the 
presence of numerous point sources of X-ray emission 
along the line of sight to MSH 11-6^. These appear to 
be associated with the open cluster Tr 18, located at a 
distance of ^ 1.5 kpc (Vazquez & Feinstein 1990). 

We have analyzed spectra from both the central X-ray 
bar and the SNR shell using background regions selected 
from regions adjacent to the bar, and just outside bound- 
aries of the remnant, respectively. Weighted response 
matrix (RMF) and ancillary response (ARE) files were 
created for the extended regions, and spectral modeling 
was then carried out using xspec version 11.3.2ag. Point 
sources within in the regions of interest were excluded 
from the spectral extraction. 

The XMM data were reduced using SAS Version 10.0. 
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Fig. 3.— Chandra ACIS X-ray spectra for PWN (black) and NS 
(red) with absorbed power law fits. 

Standard screening resulted in cleaned exposures of 35.6 
(29.2) ks for M0S2 (pn). 

3.2. Central Bar 

The central X-ray bar is aligned with the radio bar, 
and extends to the edge of the radio in the southwest, but 
the brightness falls off more quickly than the radio emis- 
sion in the northeast direction. There is a bright X-ray 
source located at the position of the brightest radio emis- 
sion, likely representing the pulsar that powers the neb- 
ula. The source, located at ll'^ll'"48".62,-60°39'26".2 
(and hereafter identified as CXOU Jill 148.6-603926), is 
roughly centered on the X-ray nebula. A fainter point 
source is detected at the southwestern edge of the radio 
bar, with coordinates ll'^ll'"44^5,-60°40'02".8 (see inset 
to Figure 2). The source is positionally coincident with 
the V = 10.8 star Tycho 8959 255 1, which is unrelated 
to the nebula. 

The Chandra spectrum of the X-ray bar (excluding 
GXOU J111148.6-603926) is well-described by a power 
law model with TV// ss 7 x lO^^ cm'^ and F « 1.7 (Ta- 
ble 1, Figure 3), securing the PWN interpretation for 
this structure. Here we have used an elliptical back- 
ground region on S3, adjacent to the PWN but within 
the SNR, to account for both sky and instrumental back- 
ground as well as thermal emission from the SNR. An 
independent fit using the XMM pn data yields virtu- 
ally identical results. The luminosity (0.5 — 10 keV) is 
1.1 X 10'^'*d| erg s""'^, about a factor of two higher than 
the value determined from ASCA observations (Harrus, 
Hughes, & Slane 1998), but far lower than the value 
of ^ 8 X lO'^^d^ erg s^^ cited by Kargaltsev & Pavlov 
(2008). 

To search for spectral variations in the PWN, we ex- 
tracted spectra from rectangular regions along the neb- 
ula, as shown in Figure 4. We fixed the column density to 
that obtained for the entire PWN, and then fit each spec- 
trum to a power law model. The results, summarized in 
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TABLE 1 

Spectral Fit Parameters 



Region 



Data 



Parameter 



Value 



PWN 

(whole) 



PWN 
Region 1 
Region 2 
Region 3 
Region 4 
Region 5 
Region 6 

NS 



SNR 
(whole) 



cxo 



cxo 



Nh 

r 

pa 



" H 



CXO 



XMM 



Ma 
" H 

r 



N°- 

net 

kT2 
^1 



(6.7 ±0.7) X 10^^ 
1.8 ±0.1 



6.7 X 10^1 

1.7 ±0.1 

1.8 ±0.1 
2.2 ±0.1 
1.5 ±0.1 
1.7 ±0.1 
2.5 ±0.4 

6.7 X 10^1 cm-2 
1.2 ±0.2 

3.7 X lO"-"^^ erg cm"^ 

6.7 X 10^1 cm-2 
l.Stoi keV 
(1.5 ±0.2) X 10^° s 
2.8 ± 0.4 keV 

2.1 X lO^-"^-"^ erg cm"^ 

3.0 X 10^-"^^ erg cm"^ 



Note. — 

a) Fixed at best-fit value for PWN 

b) Unabsorbcd 0.5 - 10.0 keV flux. 



Table 1, show a spectral steepening with distance from 
the central neutron star. Such steepening is also observed 
in G21. 5-0.9 (Slane et al. 2000), 3C 58 (Bocchino et al. 
2001; Slane et al. 2004) and other PWNe, and is un- 
derstood to be associated with synchrotron aging of the 
particles as they diffuse from the injection point to the 
outer nebula (Kennel & Coroniti 1984), although models 
matching the exact radial behavior of the spectral index 
remain elusive (e.g. Reynolds 2003). 

3.3. Compact Source 

Our Chandra observations obtain 820 background- 
subtracted counts from within a circle with a radius of 3.3 
arcsec centered on the compact source CXOU J111148.6- 
603926. The associated count rate of 1.6 x 10~^ cut s~^ 
yields negligible pileup in the image. The radial profile of 
the compact source image is well-described by the point 
spread function of the telescope, with no evidence of a 
surrounding jet/torus structure as is seen in many other 
such systems, although evidence for filamentary struc- 
ture is seen on larger scales (see Figure 4). 

The spectrum of the compact source is well-described 
by a power law with the column density fixed at the 
value obtained from fits of the entire PWN, and with F = 
1.2±0.2.The unabsorbed fiux is F^(0.5-10 keV) = 3.7x 
10"^'^ erg cm~^ s~^ indicating a luminosity ^^(0.5 — 
lOkeV) = 1.1 X 10^^ erg s~^, similar to that measured 
for a number of other young pulsars with bright PWNe 
(see Kargaltsev & Pavlov 2008). 



3.4. SNR Shell 

To investigate the X-ray properties of the SNR shell, 
we extracted spectra of the entire SNR from both the 
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Fig. 4. — Chandra image of PWN in MSH 11-62, along with 
regions used for extracting the spectral index along the nebula 
(see Table 1). The image has been smoothed with a Gaussian of 
5 arcsec width. The pulsar candidate is the bright central source 
between regions 1 and 4. White circles with red slashes indicate 
point sources removed before extraction of spectra. 

XMM pn and M0S2 detectors (see Figure 5). Emission 
from the PWN was eliminated from the spectra, as was 
emission from point sources in the field. Weighted re- 
sponse files were generated for each spectrum, and back- 
ground spectra were obtained from regions outside the 
SNR shell in each detector. 

The spectra (Figure 6) were fit simultaneously to an 
absorbed nonequilibrium ionization planar shock model 
(VPSHOCK in XSPEC), with the column density fixed at 
the best-fit value for the PWN. We find an electron tem- 
perature kT — 1.5lg'2 ionization timescale 
riet = (2.9 ± 0.9) X 10^°s cm~^, with evidence for mild 
overabundances of Ne, Mg, and Si. The fit requires a sec- 
ond, harder component, which can be fit with a power 
law with F = 2.6 ± 0.2 {xl = 1.3 with 612 degrees of 
freedom), but a slightly better fit (A^^ = 6.1 for the 
same number of degrees of freedom) is obtained with a 
second thermal component (nei model in xspec) with 
kT — 2.8 ± 0.4 keV and cosmic abundances. If associ- 
ated with the forward shock, this electron temperature 
corresponds to an expansion speed of '~ 1600 km s~^ as- 
suming temperature equilibration between electrons and 
ions (and a higher speed if equilibration has not yet been 
reached). 

Based on the volume emission measure for the lower- 
temperature component, and assuming a thin-shell mor- 
phology for the SNR, with a shell thickness of i?/12 corre- 
sponding to a shock compression ratio of 4, the estimated 

preshock density is nn = O-Oirj^^^d^^^"^ f~^^'^ cm~^, 
where 77 is the fraction of the X-ray emission that is not 
associated with point sources in the field, and / is the 
fraction of the remnant shell actually seen in X-rays. The 
observed variation in brightness around the SNR shell 
suggests possible density variations of a factor of two or 
so, and we consider the range 0.02 — 0.08 cm~'^ in mod- 
eling the dynamical evolution in Section 5. We note that 
the higher temperature component provides very nearly 
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Fig. 5.— XMM MOS2 image of MSH 11-62, with radio con- 
tours from Figure 1. The image has been exposure-corrected, and 
red/green/blue corresponds to the energy band 0.6-1.5/1.5-2.5/2.5- 
5.0 keV. The large circle indicates the region from which the SNR 
spectrum was extracted. Events from the inner ellipse were ignored 
to eliminate the contribution from the PWN, and those from the 
small circles with red slashes were ignored to reduce point source 
contributions. The same inner ellipse was used to extract the spec- 
trum of the PWN. 

the same density estimate as that from the low temper- 
ature component. 

The mass swept up by the SNR under such a sce- 
nario is then M^^ = O.brj^/^ f^/'^dl^'^ Mq. Given that the 
system morphology appears to show evidence that the 
PWN has undergone an interaction with the SNR re- 
verse shock, this low value for the swept-up mass would 
appear to suggest a rather low ejecta mass for the su- 
pernova. The inferred density, combined with the ion- 
ization timescale from the spectral fits, indicates an age 

of SSOOTy^^/^dg^^/^/^ yr, which should be considered an 
upper limit since portions of the SNR shell that are not 
detected in X-rays presumably have lower density values. 

4. GAMMA-RAY OBSERVATIONS 

4.1. Fermi LAT 

We investigated 7-ray events acquired from the region 
surrounding MSH 11-6^ with the Fermi-LAT during the 
period 2008 August 5 to 2011 March 14. Standard event 
selection was applied, using "Source" class events with 
zenith angles less than 100 degrees to minimize the por- 
tion of the Earth limb in the LAT field of view (Abdo 
et al. 2009). The Pass 7 version 6 instrument response 
functions were used, along with standard analysis tools 
available from the Fermi Science Support Center (version 
v9r21p0). 

Standard background models were used to account 
for both Galactic and extragalactic diffuse emission as 
well as instrumental background. Sources within 10° of 
MSH 11-6^, with a statistical significance larger than 
5(7 above the background, are extracted from the two- 
year Fermi-LAT Second Source Catalog (Abdo et al. 
2011), except for 2FGL J1112. 1-6040 which corresponds 
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Fig. 6.— XMM spectra from the SNR shell of MSH 11-62, with 
pn data shown in black and MOS2 data shown in red. The best-fit 
VPSHOCK model is shown as a histogram. 

to the 7-ray emission of MSH 11-6.2. The flux pa- 
rameters for all the point-like sources less than 5° from 
MSH 11-6^ are left free in the likelihood fit while the 
spectral parameters of other sources are fixed at the 
2FGL values. The mapcube file ring-2year_P76-v0.fits 
was used to describe the Galactic 7— ray emission, 
and the isotropic component was modeled using the 
isotrop-2year_P76-Source-v0.txt table. Data analysis de- 
tails follow those in Castro & Slane (2010) and Grondin 
et al (2011). 

Using the binned maximum likelihood analysis package 
gtlike, the Test Statistic, TS, (for which the significance 
is approximately {TSY^^a; Mattox et al 1996, Equation 
20) for MSH 11-6^ was found to be 432 when the source 
was analyzed in the energy range from 100 MeV to 100 
GeV. This value of the Test Statistic corresponds to a 
21 sigma detection of the source. A map of the test 
statistic as a function of position, used to localize the 
best fit position of MSH 11-6^, is shown in Figure 7. 
Overlaid on this plot are blue contours corresponding to 
ATCA observation of MSH 11-62 and white contours 
corresponding to the LAT significance levels of ~ 13, 15, 
and 17 a. 

We also indicate the positions of the three known pul- 
sars within a 45' radius from the centroid of the LAT 
source. If each of these pulsars converts its entire spin- 
down power into 7-rays in the LAT band, they could ac- 
count for no more than 6% of the observed flux, and we 
thus ignore their contributions in the subsequent anal- 
ysis. The nearest pulsar with a high E/cP ratio is 
PSR J1105-6107, located 55' from the centroid of the 
7-ray emission, well beyond the border of Figure 7; we 
consider an association with this source to be very un- 
likely. 

The source spectrum was extracted using both front 
and back events, over the energy range 0.1-100 GeV. 
The spectrum is shown in Figure 8, and is well-described 
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Fig. 7.— TS map of Fermi LAT emission from MSH 11-62. 
Blue contours correspond to the radio emission, and white con- 
tours represent the 13, 15, and 17 cr significance levels of the flux. 
Positions of nearby pulsars with E > 10^^ erg s~^ are indicated 
(see discussion in text). 

by a power law with F = 1.2 ± 0.2 accompanied by an 
exponential cutoff with Ecut = 3.2 ± 0.9 GeV that im- 
proves the fit by ~ 6.8cr. The associated flux is F{> 
100 MeV) = (6.3 ± 1.2 ± 3.2) x 10"^ photons cm"^ s'^ 
[or (8.0 ± 0.6 ± 4.2) x 10"" erg cm-^ s'^], where the 
quoted uncertainties are statistical and systematic, re- 
spectively. 

Two main systematic uncertainties can affect the LAT 
flux estimation for a point source in the Galactic plane: 
uncertainties on the Galactic diffuse background and on 
the effective area. The dominant uncertainty at low en- 
ergy comes from the Galactic diffuse emission, which we 
estimated by changing the normalization of the Galac- 
tic diffuse model artificially by 6% as done in Abdo et 
al. (2010). The second systematic is estimated by us- 
ing modified IRFs whose effective areas bracket those of 
our nominal Instrument Response Function (IRF). These 
bracketing IRFs are defined by envelopes above and be- 
low the nominal energy dependence of the effective area 
by linearly connecting differences of (10%, 5%, 20%) at 
log(E) of (2, 2.75, 4), respectively. 

The source 2FGL J1112. 1-6040, from the LAT 2"'' year 
catalog, was identified with a possible association with 
MSH ll-6^(Abdo et al. 2011). The 2FGL flux agrees 
with our measurements as well, but both are signifi- 
cantly lower than that of 3EG J1102-6103, whose quoted 
spectral index is also significantly steeper than that for 
MSH 11-6^. These differences are presumably associ- 
ated with the previously-unrecognized spectral cutoff. 

We have searched for possible extension of the source 
using the routine pointlike (Kerr 2011), comparing the 
test statistic TSext for the extended hypothesis to that 
for a point source hypothesis, TSp, defining TS^xt = 
TS — TSp. We selected events converting in the front 
section of the tracker, where the Point Spread Function 
(PSF) is the narrowest. We fitted the position of the 
source assuming iteratively a point source and a Gaussian 
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Fig. 8. — Fermi LAT spectrum of MSH 11-62. Black error bars 
correspond to statistical errors, while red error bars include sys- 
tematic errors as well. The solid curve corresponds to the best-fit 
model for a power law with an exponential cutoff. The dashed ma- 
genta curve corresponds to a pion-decay model assuming a power 
law proton spectrum with an exponential cutoff. See Section 5.2 
for discussion. 

shape for which we fitted the extension. We find that the 
best-fit centroid is located at Ilhl2ml2s, -60°40'12", 
in good agreement with the position of MSH 11-62, 
and that there is no evidence for extent with signifi- 
cance greater than Scr. Computing a significance map, 
the POINTLIKE routine determines the likelihood gradi- 
ent with spatial coordinates; this gradient has a Icr width 
of 1.2 arcmin in each direction for MSH 11-6^. The sig- 
nificance contours thus provide an indirect way of char- 
acterizing the statistical uncertainties on the position of 
the detected source. A deeper observation is required to 
firmly address, with higher confidence, whether or not 
the Fermi LAT source is extended. 

4.2. H.E.S.S. Upper Limits 

The position of MSH 11-6^ was covered in a survey of 
the inner Galactic plane with the High Energy Spectro- 
scopic System (H.E.S.S.) in 2004/2005. The approximate 
exposure time on the source was 9.1 hr (Hoppe 2008), 
resulting in no reported detection of the source. Based 
on this non-detection, we calculate a 2a upper limit on 
the VHE flux at 1 TeV for MSH 11-6^ by scaling the 
H.E.S.S. sensitivity for a 5a point source detection (1% 
of the Crab in 25 h under the assumption of a photon in- 
dex of 2.6) to the actual H.E.S.S. exposure for the source. 
We find E^{dN/dE)E=iTeV < 3.8x10-^2 erg cm'^ s-\ 
where N is the number of photons per unit area per sec- 
ond. 

5. ANALYSIS AND DISCUSSION 

5.1. Evolutionary State 

The large-scale evolution of a composite SNR depends 
on the mechanical energy of the explosion, the density of 
the ambient medium, the mass of the supernova ejecta. 
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Fig. 9. — Upper: Time evolution of the PWN magnetic field. 
Lower: Evolution of the SNR (blue) and PWN (red) radii as a func- 
tion of time. The observed radius of MSH 11-6S is indicated, and 
the vertical line indicates the age at which this radius is reached. 
See Section 5.1 for model details. 

and the spin-down power of the pulsar. Additional pa- 
rameters of considerable importance are the braking in- 
dex and spin-down timescale of the pulsar, along with the 
fraction of the spin-down power that appears as magnetic 
flux. 

Figure 9 shows an illustrative model for the evolution 
of MSH 11-62 using parameters for a scenario described 
in Section 5.3. For these parameters, the forward shock 
(FS) of the SNR (blue curve) reaches the observed radius 
of MSH 11-6^ {RsNR ~ 8.3^5 pc) at an age of - 1.2 kyr. 
By this time, the reverse shock has begun to approach the 
PWN (for which the model yields RpwN ~ 5 pc, or an 
angular diameter of « 6.9 arcmin at a distance of 5 kpc), 
as indicated by the subsequent compression of the PWN 
radius (red curve). This is consistent with the somewhat 
distorted morphology of MSH 11-6^, which appears to 
indicate that the RS has propagated most rapidly along 
the NW/SE direction. In addition, Figure 9 shows that 
the magnetic field in the PWN has declined to a value of 

The broadband emission observed from MSH 11-6^ 
places considerable constraints on the underlying proper- 
ties of the composite system. Of particular importance is 
our understanding of the observed 7-ray emission, which 
could originate from the pulsar, the PWN, or the SNR it- 
self. Here we consider models for each scenario using con- 
straints derived from the radio and X-ray observations 
of MSH 11-62. In these model investigations, we have 
varied input parameters to investigate the broad effects 
they have on the broadband spectrum, and have then 
fine-tuned parameters to provide reasonable agreement 
with the data. The resulting models are not statistical 
fits to the observed emission, and were not optimized by 
a likelihood analysis over the large parameter space. 

5.2. Gamma-rays from the SNR 



Particle acceleration at the forward shocks of SNRs 
is known to produce extremely energetic particles that 
can produce 7-rays through inverse-Compton (IC) scat- 
tering of ambient photons by energetic electrons, and/or 
through the decay of neutral pions produced in energetic 
proton-proton collisions. Thus, one interpretation for the 
7-ray emission from MSH 11-62 is that it originates from 
the SNR itself, and is a signature of particle acceleration 
by the remnant. We consider two cases, one in which the 
7-ray emission is dominated by hadrons (vr'^-decay) , and 
one in which it is dominated by leptons (IC scattering). 

For the hadronic-dominated scenario, we assume that 
the mechanical energy of the SNR is Esnr = 10^^ erg 
and that 50% of this energy is converted to energetic par- 
ticles whose spectrum is described by a power law with 
an exponential cutoff. We assume an electron-to-proton 
ratio ke-p = 10~^, and that the spectral index is the 
same for protons and electrons. Using d = 5 kpc, we 
adjust the spectral index, density, and cutoff energy for 
the proton spectrum to match the Fermi LAT spectrum, 
while also adjusting the magnetic field strength to repro- 
duce the SNR radio emission with synchrotron emission 
from the electron spectrum. We find a spectral index 
of 1.8 and an exponential cutoff energy of 70 GeV for 
the protons (Figure 10, top), with an ambient density 
no = 6.8 cm^'^. The associated 7-ray emission from IC 
scattering of the cosmic microwave background (CMB) 
is negligible. The postshock magnetic field required by 
the radio emission is — 10/iG. Here we have used the 
measured 1.4 GHz flux from the shell, and assumed a 
spectral index range of = 0.4 — 0.7 typical of radio 
emission from shell-type SNRs. 

The required proton cutoff energy under this scenario 
is rather low given that we expect energies well in ex- 
cess of 10^^ eV to be produced in young SNRs, but it 
is possible that a low shock speed has resulted in a low 
maximum energy (e.g., Reynolds 2008). Of potentially 
greater concern is that the preshock density required un- 
der this scenario is more than two orders of magnitude 
larger than that inferred from the X-ray spectrum. Sim- 
ilar discrepancies have been identified for several 7-ray 
emitting SNRs associated with OH masers (which im- 
ply an interaction between the SNRs and surrounding 
molecular clouds, thus making Tr'^-decay a likely produc- 
tion mechanism), and it has been suggested that this 
could be the result of a clumpy SNR shell in which only 
the lower-density interclump medium is responsible for 
the X-ray emission while all of the gas participates in the 
7-ray production (Castro & Slane 2010). No maser emis- 
sion has been reported for MSH 11-6^, and the remnant 
is not known to be interacting with dense material. Thus, 
while this scenario cannot be ruled out, particularly since 
the production of masers requires a rather narrow range 
of physical conditions, it does not have strong supporting 
evidence outside of the 7-ray emission. Searches for evi- 
dence of interactions between MSH 11-62 and molecular 
clouds are of particular importance to further address 
this picture. 

For at least some SNRs, the observed 7-ray emission 
appears to be dominated by IC emission from an ener- 
getic population of electrons (e.g., Elhson et al. 2010). 
We have considered such a scenario for MSH 11-62 by 
fixing the ambient density at the value implied by the X- 
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Fig. 10. — Models for the broadband emission from MSH 11- 
62 under the assumption of dominant hadronic (upper) and lep- 
tonic (lower) processes for the observed 7-ray emission. The dotted 
curves correspond to synchrotron emission, and the dashed curves 
correspond to IC emission. The red curve represents 7-ray emission 
from pion decay, and the solid black line represents the summed 
7-ray components. 

ray measurements, and allowing the normalization and 
shape of the electron spectrum to vary. We find that for 
a distance of 5 kpc, the total electron energy required 
to produce the 7-ray emission through IC scattering of 
CMB photons is 8 x 10^^ erg. The magnetic field re- 
quired to fit the radio emission from the SNR shell is 
B = 2.0/LtG. The model, shown in Figure 10 (bottom), 
requires a super-exponential cutoff of the electron spec- 
trum 



Ne oc E-°'' exp 



En 



(1) 



Model parameters are given in Table 2. In particular, we 
find Ecut.e = 0.9 TeV, which is quite low given the rather 
low magnetic field. Since all realistic acceleration models 
that produce energetic nonthermal electrons in SNRs also 
accelerate ions, the total energy budget implied by these 
results would exceed the available energy of the SNR 
blast wave for fcg-p ^ 10~^ if the remnant distance is 
5 kpc. 

For the smallest distance allowed by HI absorption 
measurements, an electron-to-proton ratio fce-p > 0.7 
is required for the entire particle energy budget to not 
exceed 5 x 10^^ erg. Such high required values of fce-p, 
along with the low cutoff energy and extremely sharp 



TABLE 2 

Model Parameters: 7-rays from SNR 



Parameter 


Hadronic Model 


Leptonic Model 


d (fixed) 


5 kpc 


5 kpc 


no 


6.8 cm-3 


0.04 cm-3 


OLe {= ap) 


1.8 


1.8 




50 GeV 


0.9 TeV 




70 GeV 


1 TeV 




1.0 


2.0 


El 


5 X 10** erg 


8 X 10*^ erg 




5 X 10^° erg 


8 X 10^5 1 erg 


fce-p (fixed) 


10-2 


10-2 


B 


lO.OAtG 


2.0aiG 



Note. — 

a) Curvature index on exponential eut-off 

b) Total electron energy 

c) Total proton energy 

spectral cutoff appear problematic. Thus, we conclude 
that a leptonic scenario in which the observed 7-ray emis- 
sion from MSH 11-6^ is dominated by IC emission is 
ruled out for all but the smallest possible distances, and 
still problematic at such small distance values. 

5.3. Gamma-rays from the PWN 

The X-ray and radio observations of MSH 11-62 
clearly identify a powerful PWN in the system. The elec- 
trons responsible for the X-ray and radio emission also 
produce 7-ray emission through inverse Compton scat- 
tering. To investigate the scenario in which the Fermi 
LAT emission is produced entirely by the PWN, we cal- 
culate the evolution of the composite system using the 
model of Gelfand, Slane, and Zhang (2009). We assume 
expansion of the SNR into a uniform medium of density 
riQ, driven by an ejecta mass Mej. The PWN expan- 
sion is powered by the input of a central pulsar with a 
spin-down power that evolves as 



E{t) - Eo 



1 



t 

To 



P+1 
p-1 



(2) 



where tq is the initial spin-down timescale for the pulsar, 
and p is the braking index. The PWN expansion is con- 
fined by the surrounding ejecta. Particles are injected 
into the nebula with a power law spectrum, and evolved 
through radiative and adiabatic losses. The emission 
spectrum is then calculated based upon the evolved mag- 
netic field strength of the nebula at the age required to 
reproduce the observed radius of the SNR (Figure 9). 
The ambient radiation field was assumed to be a combi- 
nation of the cosmic microwave background (CMB) and 
infrared photons from local dust. 

The parameters of the model are summarized in Ta- 
ble 3. We fixed the distance at 5 kpc and assumed 
EsNR = 10^^ erg. We sampled values for tiq in the range 
0.02 — 0.08 cm-^ based on the estimates provided from 
the X-ray emission (see Section 3.4), and determined that 
the broadband spectrum is not highly sensitive to the 
ambient density, with fiux values varying by only ~ 30% 
in the X-ray band (see Figure 11). Braking index values 
over the range p = 2 — 3, and ejecta mass values between 
1.5 — 8.5Mq, produce flux variations similar in magni- 
tude to those observed for the above range in no, and we 
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TABLE 3 

Model Parameters: 7-rays from PWN or Pulsar 



10 11 12 13 14 15 

Log Electron Energy, E, (eV) 
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Fig. 11. — Models for the nonthermal emission from MSH 11— 
6S (lower panel) under the assumption that the 7-ray emission is 
produced entirely by the PWN. Dashed, solid, and dotted lines 
correspond to no = 0.02, 0.04, and 0.08 cm~^. The upper panel 
shows the associated electron distribution, which is modeled as an 
evolved Maxwellian distribution with power law tail. 

thus fixed these quantities at the values shown in Table 3. 
The observed radio flux provides an upper limit for the 
minimum electron energy; lower values of Ee^min do not 
appreciably change the spectrum. Similarly, the X-ray 
spectrum provides a lower limit for the maximum elec- 
tron energy. Higher values of E^^max allow the photon 
spectrum to extend to higher energies, but the current 
data provide no additional constraints on this. 

With the above parameters fixed, we varied values for 
the pulsar spin-down timescale, tq, the initial spin-down 
energy, Eq, the electron spectral index, ag, and the frac- 
tion of the pulsar wind energy in the form of magnetic 
fields and particles, rjs and r/e respectively. We found 
that it was not possible to reproduce the significant 7-ray 
emission in the Fermi LAT band with simple power law or 
broken power law models for the particle spectrum under 
these assumptions. We then investigated a model con- 
sisting of a Maxwellian electron distribution with a power 
law tail, corresponding to spectra observed in PIC sim- 
ulations of relativistic shocks (Spitkovsky 2008), consid- 
ering a range of values for the energy of the Maxwellian 
peak and for the fraction of the total electron energy 
found in the power law component. The temperature 
and energy density of the IR photon field were adjusted 
slightly to tune the agreement with the data. The final 
values listed along with the other model parameters in 
Table 3 are within reasonable expectations for this region 
of the Galaxy (see Strong et al. 2000). 

As shown in Figure 11, we find that the broadband 
emission from MSH 11-6.2 can be reproduced to reason- 
able fidelity for the parameters shown in Table 3. The 
upper panel of the Figure contains the electron spec- 
trum, and the lower panel presents the synchrotron and 
IC emission along with the observed data. The radii for 
the SNR and PWN are shown in Figure 9, along with 
the PWN magnetic field. The age required to match 
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Input: 






d (fixed) 


5 kpc 


5 kpc 


EsNR (fixed) 


10^1 erg 


10^1 erg 


(fixed) 


2.5Mq 


2.5M0 


no (fixed) 


0.04 cm-3 


0.04 cm-3 


p (fixed) 


3.0 


3.0 


TO 


1000 yr 


1000 yr (fixed) 


Eo 


4 X 10^^ erg s~^ 


3 X 10^** erg s"! 




0.999 


0.995 


'IB 


0.001 


0.005 


TP 


0.1 GeV 


1 MeV 




900 TeV 


8.5 PeV 




100 GeV 




J PL 


.1% 




J^l 




4.4 GeV 


OLe 


2.33 


0.0 


ae,2 




2.7 


rpi 


25K 


25K (fixed) 


(Wir/Wcmb)'' 


2WcMB 


2WcMB (fixed) 


Derived: 






Age 


1200 yr 


1300 yr 




e.lfiG 


7.3AtG 



Note. — 

a) Fraction of energy in electrons 

b) Fraction of energy in magnetic field 

c) Peak energy of Maxwellian electron component 

d) Fraction of energy in power law component 

e) Power law break energy 

f) Effective temperature of IR photon field 

g) Ratio of IR energy density relative to CMB 



the observed SNR radius is ~ 1.2 kyr. The magnetic 
field of the PWN has declined to a value of ~ 6/iG, and 
the PWN radius is approaching the stage at which it 
encounters the SNR reverse shock, as indicated by the 
subsequent decline in RpwN- This appears qualitatively 
consistent with the X-ray and radio morphology, as dis- 
cussed in Section 5.1, however the narrowest portion of 
the PWN is considerably smaller than indicated in Fig- 
ure 9. The very distorted shape of the nebula suggests 
that the reverse shock has propagated asymmetrically in 
MSH 11-6^, which is not a scenario that can be treated 
in our spherical model. 

The magnetic field in this scenario is far below the 
minimum-energy value calculated from radio measure- 
ments, indicating that the nebula is particle-dominated. 
We note that the ^ 65/zG field estimated by Harrus, 
Hughes, & Slane (1998) was based on the assumption 
that the minimum-energy conditions hold in MSH 11- 
62. Additional estimates based on the interpretation of 
the spectral break between the radio and X-ray bands as 
being an evolutionary effect, and on the pressure in the 
PWN being in equilibrium with that of the thermal gas 
in the SNR, led to even higher estimated values of the 
magnetic field. However, no modeling was performed to 
demonstrate which, if any, of those assumed conditions 
provide a self-consistent picture of the composite sys- 
tem. Here we have carried out such modeling efforts and 
find that a particle-dominated nebula is required, with 
the magnetic field at the current epoch being well be- 
low the equipartition value. Similar results are obtained 
for the evolved nebulae in W44 (Bucciantini et al. 2011), 
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GO.9+0.1, and G338.3-0.0 (Fang & Zhang 2010). We also 
note that even though different values for the assumed 
distance can lead to different values of the model pa- 
rameters (see below), the observed X-ray spectral index 
constrains the magnetic field to be low, and the overall 
magnetization to be low as well. 

According to this model, 7-ray emission in the Fermi 
LAT band is produced entirely by electrons in the 
Maxwellian peak of the spectrum. This is consistent 
with results for Vela X (de Jager, Slane, & LaMassa 
2008), HESS J1640-465 (Slane et al. 2010), and sev- 
eral other PWNe (Fang & Zhang 2010) that suggest a 
low-energy electron component distinct from the power 
law component (although for Vela X this was modeled as 
an additional power law component). For MSH 11-62, 
we find a mean Lorentz factor of 7 = 2 x 10^ for the 
electrons in the Maxwellian component, with only 0.1% 
of the total energy in the power law tail. 

We note that the initial spin-down power required for 
this model, Eq = 4x 10^^ erg s~^, is uncomfortably large. 
For the assumed braking index p = 3, this would lead to 
a value of £^ = 8.3 x 10^^ erg s~^ at the current epoch 
- higher than that for any known pulsar. This would 
imply Lx/ E ^ 10^^ for the pulsar, which is extremely 
low. 

Given that the distance to MSH 11-6^ is not well 
known, the numerical values from the model above are 
quite uncertain. Changes in the assumed distance can be 
accommodated, at least in part, by adjusting the spin- 
down power, for example. However, larger distance val- 
ues lead to even higher values of E. Reasonable fits 
to the broadband spectrum can also be obtained for a 
distance as low as 2 kpc, yielding a current spin-down 
power of 2 X 10'^^ erg s^^. While still extremely high, 
this is similar to the value for PSR J20224-3842, a pul- 
sar recently discovered in G76. 9-1-1.0 (Arzoumanian et 
al. 2011). This pulsar is characterized by an extremely 
low conversion efficiency of spin-down power into X-rays 
as well, but the value {Lx/E ~ 6 x 10^'°^) is still much 
higher than that inferred for the pulsar in MSH 11-62 if 
the Fermi LAT emission is being produced by the PWN. 
We conclude that it is possible for the observed 7-ray 
emission from MSH 11-6^ to be produced by the PWN, 
but that this would require a very significant amount of 
energy in the Maxwellian peak of the electron spectrum, 
leading to spin-down parameters for the pulsar that seem 
unlikely. 

5.4. Gamma-rays from the Pulsar 

Fermi observations have clearly established that pul- 
sars comprise a significant fraction of the Galactic popu- 
lation of GeV point sources. The associated spectra are 
characterized by a power law with an exponential cutoff, 
with E'cut -^1-6 GeV, and F ~ 1 - 2 (Abdo et al. 
2010). This is very similar to the spectrum observed for 
MSH 11-62. Being associated with a composite SNR, 
for which our X-ray observations identify the underlying 
neutron star, it is quite possible that the observed 7-ray 
emission originates with the pulsar itself. 

The ratio of the X-ray fiux from the point source in 
MSH 11-6^ to the observed 7-ray flux is F^/F^ 5 x 
10~^. This value is typical of that for pulsars detected 
by the Fermi LAT, although there is a large spread in 
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Fig. 12.— Model for the PWN emission from MSH 11-6S un- 
der the assumption that the 7-ray emission arises from the pulsar. 
Parameters for the model are given in Table 3. 

this ratio for the population. 

Using the approximate relationship ~ IQ^^E de- 
rived from a variety of studies of X-ray PWNe and 
their associated pulsars, the spectrum of the PWN in 
MSH 11-6^ suggests an underlying pulsar with E ~ 
lO'^^d^ erg s~^. We note, however, that there are large 
variations in the observed relation between and E 
(e.g., Vink, Bamba, & Yamazaki 2011). More than 85% 
of the pulsars detected by the Fermi LAT, and all of 
those pulsars discovered in the 7-ray band, have lower 
values of E. However, of currently reported pulsars dis- 
covered in 7-rays, only one has a Fermi LAT flux lower 
than that associated with MSH 11-62. Thus, the lack of 
identified pulsations in any band is not strong evidence 
that the emission from MSH 11-6^ is not associated with 
the pulsar itself. The marginal evidence for extension of 
the 7-ray emission, which would rule out the pulsar as 
the source of the emission, is of particular importance in 
this context, and requires revisiting in the future when 
additional data from this source have been obtained. 

If we assume that the Fermi LAT emission arises en- 
tirely from the pulsar in MSH 11-6^, the PWN radio 
and X-ray emission can then be reproduced by a broken 
power law electron spectrum with i?o 3 x 10'^* erg 
and an age of 1300 yr, at which time the ambient 
density tiq — 0.04 cm~^ inferred from the SNR X-ray 
emission yields the observed radius for the SNR assum- 
ing a distance of 5 kpc. Here, as in Section 5.3, we have 
fixed the values for rf, M^j, no, and p. In addition, we 
have used the same model values for the IC photon fields 
as were used in Section 5.3, although these are largely 
irrelevant here since the IC emission falls well below the 
observed 7-ray emission. We have varied the remaining 
parameters in Table 3 in order to provide an approx- 
imate match to the observed radio and X-ray spectra. 
The resulting model spectrum is shown in Figure 12, 
and the parameter values for the model are listed in Ta- 
ble 3. The PWN magnetic field strength is ~ 7/xG, a 
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factor of 5 lower than the minimum-energy field derived 
from radio observations, and the spin-down power at the 
current epoch is ii^ = 5.7 x 10'^'' erg s~^. This leads to 
Lx/E « 2 X 10~^, somewhat low, but well within the ob- 
served range for known pulsars. In this model, the SNR 
RS has not yet begun to impact the PWN, leaving the 
observed X-ray and radio morphology unexplained. 

We note that this model again leads to a particle- 
dominated nebula, with the magnetic energy comprising 
only ^ 4% of the total energy. Models that produce 
larger magnetic fields at the derived age can be obtained 
by increasing rjB- However this leads to a synchrotron- 
loss cutoff below the X-ray band. Conversely, a better fit 
to the X-ray slope can be obtained by further reducing 
the fraction of wind energy that appears as magnetic flux. 
This, however, requires a much higher value of Eq and 
results in an even lower value for the current magnetic 
field strength in the PWN. In any event, the low magne- 
tization of the nebula appears robust in these models, al- 
though we caution that our spherically-symmetric semi- 
analytical calculations do not address more complex geo- 
metrical and magnetohydrodynamical effects that surely 
complicate the magnetic field structure in PWNe. 

6. CONCLUSIONS 

Radio observations of MSH 11-6^ indicate that this is 
a composite SNR containing a central PWN. X-ray ob- 
servations confirm this, with clear identifications of the 
PWN synchrotron emission as well as the presence of 
an X-ray point source that most likely corresponds to 
the pulsar driving the system. The X-ray spectral in- 
dex in the PWN increases with distance from the com- 
pact source, as expected from synchrotron losses, and 
the strongly asymmetric morphology of the PWN sug- 
gests an interaction with the SNR RS. 

Observations with the Fermi LAT reveal distinct 7- 
ray emission from MSH 11-6^, but the limited angu- 
lar resolution makes it impossible to ascertain the ori- 
gin of the emission directly. We have considered three 
cases, with the 7-ray emission arising from the SNR, the 
PWN, and the pulsar, respectively. Models for each sce- 
nario can reproduce basic features of the observed broad- 
band nonthermal emission, but not without difficulties. 
If the emission arises from the SNR, then a hadronic 
scenario requires an ambient density that is a factor of 
^ 100 higher than that implied by the thermal X-ray 
emission, while a leptonic scenario requires that an un- 
reasonably large fraction of the SNR mechanical energy 
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